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Resumo 
Este trabalho apresenta uma aproximação da regra da cadeia da 
derivada fracionária, do tipo Caputo com baixo nível fracionário, 
para obter soluções exatas das versões fracionárias das equações 
diferenciais de Bernoulli e de Riccati. 
Palavras-chave: Cálculo Fracionário e Aplicações. Bernoulli. 
Riccati. Regra da Cadeia. Baixo Nível Fracionário. 


Abstract 
This work presents an approximation of the Fractional derivative 
chain rule of the Caputo type with low fractional level to obtain 
exact solutions of the fractional versions of the Bernoulli and Ric- 
cati differential equations. 
Keywords: Fractional Calculus and Applications. Bernoulli. 
Riccati. Chain Rule. Low-Level Fractionality. 
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1 Introdução 


O cálculo de ordem não inteira, conhecido como cálculo fracionário, tradução livre de fractio- 
nal calculus, generaliza o cálculo integral e diferencial clássicos. Pode-se pensar nos operadores 
de ordem fracionários como os operadores que representam funções da memória sobre a história 
de alguns sinais de sistemas físicos (CAMARGO; OLIVEIRA, 2015; PODLUBNY, 1999). 

Nos últimos anos, a caracterização de sistemas reais por meio de modelos com equações dife- 
renciais de ordem fracionária vem aumentando consideravelmente, pois os mesmos representam 
com maior fidelidade fenômenos naturais que os modelos de ordem inteira. Mais ainda, em 
relação ao tempo esses modelos descrevem processos com memória de longo prazo do tipo lei de 
potência (DU; WANG; HU, 2013; EDELMAN, 2014; JUCHEM NETO; FONSECA, 2018; PO- 
DLUBNY, 1999; TARASOV, 2013; TARASOVA; TARASOV, 2017). Desta forma, a evolução 
do sistema estudado depende de toda a memória dos eventos passados. 

Recentemente, foi proposto em Tarasov e Zaslavsky (2006), e de forma independente em 
Varalta, Gomes e Camargo (2014), casos em que a se desvia minimamente de um valor in- 
teiro wu =n-—e. Para0< e << 1 o comportamento da derivada fracionária é denominado de 
baixo nível fracionário. Essencialmente, o caso baixo nível fracionário é tradado como uma 
perturbação para o caso inteiro através da £-expansão por Tarasov e Zaslavsky (2006). Este 
método também foi aplicado a oscilação fracionária (TOFIGHI, 2013), fenômenos de relaxa- 
mento fracionário (TOFIGHI; GOLESTANI, 2008) e física das partículas (GOLDFAIN, 2008). 

Utilizamos o caso de baixo nível fracionário nas versões fracionárias das equações diferen- 
ciais dos tipos Bernoulli e de Riccati. Neste caso, podemos aproximar a derivada fracionária de 
Caputo e utilizar uma versão aproximada da regra da cadeia para obter uma solução semelhante 
a solução clássica. 

Parte deste trabalho foi apresentado durante o Encontro Regional de Matemática Aplicada e 
Computacional que ocorreu de seis a oito de junho de 2018 na Universidade Estadual Paulista 
“Júlio de Mesquita Filho” no campus de Bauru. 


2 Preliminares 


Nessa seção introduzimos a definição da derivada fracionária segundo Caputo (1967) e algu- 
mas de suas principais propriedades. Maiores detalhes podem ser vistos em Camargo e Oliveira 
(2015), Diethelm (2010) e Silva, Moreira e Moret (2018b). 











Definição 2.1 Sejam T, «e Rten=minfke N/k > a). A derivada fracionária de ordem a 
no sentido de Caputo (à esquerda) da função f E AC"[0,T] é difinida como 





t 
o os no n—1l—-a (n) d 1 Ee 
SDf(t)=4 T(n-o) A (x —s) f”(sjds, n n, 
fe), dE 
onde ft (s) = d" f(s)/ds”" eI(:) éa função Gamma. Usaremos a notação simplificada Dº f(t) = 
COCrE), 
0% 


Algumas propriedades da derivada de Caputo são bastante semelhantes às da derivada clássica 
(Newton-Leibniz). A seguir apresentaremos algumas propriedades: 


(1) 


Propriedade 1 D?K =0, para toda constante K. 
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Propriedade 2 DHAf(t)+g(t)]=A DP Ff(t) + DE g(t), onde À ER e f,g E AC"[0,7]. 
Teorema 2.2 (DIETHELM, 2010) Seja f(t) = t*, com À >0, então temos a seguinte relação 


arn- DA+D qa 


Teorema 2.3 (LI; DENG, 2007) Sejam fe C""!o,T]ea =n-—e, então 
(a) lim DIE f(1) = 100 (1), 
E507 
(b) tim DEEF(O) = f0(1) - ft (0). 
E>0* 
No entanto, existem algumas inconsistência tais como: 


(1) Não obedece à familiar Regra do Produto para duas funções: 
DFL (O) - el) £ DP FCO) (1) + fl) - DP g(t). 
(2) Não obedeçe à Regra de Quociente para duas funções: 


aTf(DT , e(DDEF(O)) — f(x) DE (e(t)) 
Pr | | 7 gê(t) 





(3) Não obedece à Regra da Cadeia 
DFL f(et)] £ Dãaf(e(t))- DP e(t). 


Segundo Tarasov (2016) essas particularidades nos permitem descrever novos tipos incomus 
de sistemas físicos, naturais ou não. A violação da regra da cadeia é uma das principais carac- 
terísticas das derivadas fracionárias (SILVA, 2018). 

Deve-se enfatizar que existe uma complicada regra da cadeia para o operador de Caputo, a 
Fórmula de Faà di Bruno (DIETHELM, 2010), sendo considerada uma limitação prática deste 
operador, e de outros operadores de ordem fracionária (TARASOV, 2016). 

Neste trabalho, consideramos a seguinte equação diferencial fracionária envolvendo a deri- 
vada de Caputo: 

Dt) =f(ta(t)), a E (0,1), (3) 


com a condição inicial xo = x(0). 


Definição 2.4 A função de Mittag-Leffler de um parâmetro é dada por 
H 


=) = (4) 
= oa (ak+ 1) 
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Assim como a natureza da função exponencial é a solução de equações diferencias de ordem 
inteira, a função Mittag-Leffler tem um papel análogo para solução de muitas equações diferen- 
ciais de ordem não inteira. No início do século passado, Gôsta Magnus Mittag-Leffler apresentou 
a função inteira 


eo n 


PA 
Ea(z) = E Tok) ze C,R(z) >0, (5) 


que leva seu nome, e investigou suas propriedades básicas (MITTAG-LEFFLER, 1903, 1904, 
1905). Para 0 < a < 1, ela interpola a função exponencial e a função hipergeométrica 15 E 
fácil ver que a função de Mittag-Leffler é uma generalização da função exponencial, pois quando 


o = 1 temos a E 
El) =) —— = exp(2). (6) 
É x P(k+1) 




















Além disso, para a E (0,1), a restrição de Eg : R — R é estritamente crescente, como obser- 
vado em Cong, Son e Tuan (2014), Hilfer e Seybold (2006) e Pollard (1948) com 
dim Eo ()=o e dim Ea (1) =0. (7) 
A função de Mittag-Leffler foi generalizada por Wiman (1905), que apresentou a função de 
Mittag-Leffler com dois parâmetros, a e B, da forma 


Ze 


Essl)=L ak rp) (8) 
com a,B EC, R(a) > 0. Também é conhecida como função de Wiman. Estas funções permane- 
ceram quase esquecidas durante mais de cinguentas anos, até que foram retomadas por Agarwal 
(1953), Humbert (1953), e vários outros. 

Como caso especial, a função de Mittag-Leffler com um parâmetro é obtida para B = 1, isto 
é, Ea 1(z) = Ea(z). Outra relação importante entre as funções Mittag-Leffler com um e dois 
parâmetros é dada por (PODLUBNY, 1999) 


[Eo(z) — 1). (9) 


Eg,a+1 (z) EAST 
Pé 


2.1 Transformada de Laplace 


A transformada de Laplace é um conceito relevante para a resolução de equações envolvendo 
derivadas e integrais fracionárias. A seguir estão listadas algumas propriedades importantes para 
este trabalho. As demonstrações podem ser encontradas em Camargo e Oliveira (2015), Silva, 
Moreira e Moret (2018a): 





LIpfr(m]=s"F(s)-s"!'H(0), 0<a<i1, (10) 
s sa 
L Eo(-Atº)] A (1) 
É B 
E PE 4a 
SIA =1"" Ea pl- MO). (12) 
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3 Regra da cadeia de ordem « = 1 - e, come < 1 


Considerando n = 1 e « = 1 — £ na Equação (1) obtemos 











DI-tf(g) = = / Sia (13) 
= S. 
É F(e) Jo (s—t)l=e 
Assumindo f € C2[0,T] e aplicando o Teorema (2.3) para t E (0,T), obtemos 
' — Anflt) 
lim DIE F(t) = f'(t) = lim 14 
dim DE CH) = 0) = lim =, (14) 
onde As f(t) = f(t+e) — f(t) denota o operador diferença progressiva. 
Para a próxima definição precisaremos da seguinte equação lim, ,0+ Ane) = lim, 0 A fto) = 
lim, s0- Saito conhecida do cálculo clássico. 
Assim, podemos escrever a equação (14) na forma 
= o Meo) Acf) 
lim D/*H(N)=1 Eos o(£), 15 
er Ê H Pe € € E Co) 


em que O(£) é o erro de truncamento de ordem £ (£ > 0). 
Dessa forma, a aproximação de primeira ordem para a derivada de Caputo com baixa fracta- 
lidade,0<e< 1: 
Acf(t) 


l-e Ay 
ERC Io E (16) 





com erro de ordem e. 
As Fig. (1) e Fig. (2) exibem o gráfico das funções f(t) = DI “(t*) e g(t) = DI “(eN), 


PG) 3+e e 24º E 1+e(2t), respecti- 


utilizando Eqg. (16), para e = 0.01. As soluções exatas são r(a+e) 


vamente. 
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Figura 1: Solução exata versus solução aproximada da função 22º? (1º). 


Exata 
— — — Aproximada 








Esta aproximação de primeira ordem é válida para todos os operadores fracionários que pos- 
suem a propriedade da aproximação dada pelo Teorema (2.3). 

Quando o incremento £ é pequeno o suficiente, mas constante, podemos dizer que a diferença 
finita progressiva pelo incremento tende a derivada, mas não é igual a derivada de Caputo. Da 
mesma forma que o lado direito da Equação (16) também não é a derivada clássica. 

Aplicando o Teorema (2.3) na função composta f og, com fegemC? [0,7], tem-se 


Jim, Di [o g)(0) = (De fl) eg 80)» (17) 
 [Acf(g) |. | Aeg(t) 
= lim, to a [E (8) 


Imediatamente obtemos três aproximações possíveis combinando D/"* ou D! na Equação (18) 
com £ < 1: 


Di fog) = (Dele) agr Cel), (19) 

Di lfog(t) = (Dj es) «Digo (20) 
e 

Diogo) = (DIE) mei CE (O). (21) 


Escolhendo a Equação (19) como derivada de Caputo com baixo nível fracionário, o — 1, nós 
obtemos 


DE (Fe(1))) = (Def (8) eq) DES (O) - (22) 
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A regra da cadeia, derivada de ordem inteira, para a derivada de funções compostas dá origem 
ao chamado método da substituição ou de mudança de variável para resolver equações diferen- 
ciais. Esta técnica requer pouco cálculo e converte algumas equações diferenciais em tipos que 
podem ter suas soluções exatas conhecidas. Por exemplo, renomeando uma variável dependente 
para que u(t) = y(z(t)). Segue-se que 


Diu=y (2) Dya. (23) 


Com a escolha adequada de z, em situações especiais, a equação transformada (23) pode ser uma 
equação que possua solução exata conhecida. 


4 Resultados 


Usando a nova regra da cadeia para derivada de Caputo com baixo nível fracionário estudamos 
a equação diferencial fracionária dos tipos Bernoulli e o de Riccati. 
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4.1 Equações diferenciais fracionárias de Bernoulli com baixo nível fra- 
cionário 
Consideremos a versão de baixo nível fracionário para a clássica equação diferencial de Ber- 
noulli com O < £ < 1, dada por: 


À Dj y(t) + B(t)y(t) 
v(0) 


onde A(t) e B(t) são funções contínuas e q uma constante. Deve-se notar que algumas condições 
suficientes para a existência global de soluções da Equação (24) são dadas por Báleanu e Mustafa 
(2010) e Lakshmikantham, Leela e Devi (2009). 

Para este fim, alteremos a variável dependente de y para u usando a seguinte substituição 


u(t) = [(o)]'*. (25) 


A(Dyi(t) 
Jo Es 0, Ee) 


Consequentemente, em vista de (19) temos, 
D; “ult) = (1 q)b(t)] ID, Fu(t). (26) 
Observe que se q £ 0 e q 1, a Equação (24) pode ser reduzida a uma forma linearizada 
Dy Fut) + (1 q)B(u(t) = (1-q)A(r), (27) 


cuja solução exata pode ser encontrada, em alguns casos, usando transformada de Laplace (ATA- 
NACKOVIC et al., 2014; CONG; SON; TUAN, 2014; LIANG; WU; CHEN, 2015). 
Se para cada índice £, us; é uma solução da Equação (24) e se o limite 


y(1) = lim [ue(t)] Pá (28) 


e50+ 


existe para todo t E [0,T], então y é solução do problema no caso £ = 0, isto é, para equação de 
diferencial ordinária de Bernoulli. 

Para encerrar esta seção apresentaremos dois exemplos. O último apresenta uma função 
exponencial deformada com dois parâmetros. 


4.1.1 Equações diferenciais fracionárias de Bertallanfy-Richards com baixo nível fracionário 
Considere a Eg. (24) com A(t) = a, B(t) = b e a condição inicial y(0) = yo > 0. 
Di xl) =ay(t)-by(t), 0<e<a, (29) 


onde a 0 e b £ O são constantes arbitrárias. Ao aplicar a mudança de variável conforme definido 
em (25), a Eg. (29) é tranformada em 


D'*u=-(1I-gbu+(I-ga. (30) 
Aplicando a transformação de Laplace para ambos os lados, segue 


shºU(s)=su(0)-(1-g)bU(s)+(1-g)bs |. (31) 
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Ou seja, 
()=(1-gaz” e +u(0) 2" (32) 
ida da sie +(1—g)b sie +(1-g)b] 
A partir disso, pela Eg. (12), temos: 
u(t) = (1 g)at Er e 2-el(g — bt] + uoEr-ellg — 1)br' "8. (33) 


Então, usando a Eg. (9) e Eg. (25), de onde conclui-se 


v(t) = [5 + (uy! 5) Ei-ella— Dbr'=*]] rã (34) 


Se (q — 1)b < 0, da Eq. (7), a solução (34) converge para 


p= vie) = (8) (35) 


Contudo, como indicado em Varalta, Gomes e Camargo (2014) se € — 0 então a função dada 
na Eg. (29) é equivalente a equação diferencial de Bertallanfy-Richards. 


v(t) = | + Cm = -) eba RA (36) 


Exemplo 1 Em particular, para discutir o relaxamento não-linear Tofighi (2012) considerou 
a=1,b=2ey(0) = 1 na equação (29), 


Di*y4+2y-y =0. (37) 


A solução da Eq. (37) com baixo nível fracionário é dada por 





2 
v(t) — 1+Ey-e(2t1-8) . (38) 
4.1.2 (q,0)-Exponencial 
Considere Eq. (24) com A(t) =a £ 0, B(t) = 0 e a condição inicial y(0) = 1. 
DP vt) = ay“(t), (39) 


onde a é uma constante arbitráriae «= | —- e, com0< e < 1. Ao aplicar a mudança de variável 
conforme definido em (25), a Eg. (39) é tranformada em 


Di“u=(1I-ga. (40) 


Análogo a solução da Eq. (29), segue que 


(1) = [! tarde e (41) 
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4.2 Equações diferenciais fracionárias de Riccati com baixo nível fracionário 


As equações diferenciais com baixo nível fracionário e do tipo Riccati são da forma 


Dip dB +c(t) = A(Dy? 
À v(t) + Blt)y(t) ao E o (42) 


Se C(t) = 0, a equação diferencial de Riccati reduz-se a um caso particular da equação diferen- 
cial fracionária de Bernoulli com q = 2. No caso de A(t) = 0, recai numa equação diferencial 
fracionária do tipo linear. 

Se A(t) £ O e conseguirmos de alguma forma obter uma solução particular da Eg. (42) yp, 
então vale o seguinte resultado. 


Teorema 4.1 Seja yp uma solução particular da Equação (42). Então a substituição 
= 
v(8) =yp(t) + [z(t)] (43) 


transforma a equação diferencial fracionária de Riccati num equação diferencial linear fra- 
cionária na variável z 


Dj *z(t) + [2A(t)yp(t) — B(B)Je(t) = —A(1). (44) 


Demonstração: Da Eq. (43) obtém-se 








l-g 
Dj *y(t) = Di *yp(t)— o (45) 
Substituindo em (42), obtemos 
= pr B 2A A 
DES + Bolo) +) ET = (ob (PR, as. (49) 


Como yp é uma solução particular de (42), D/ *yp(t) + B(t)yp(t) + C(t) = A() [yp(t)]? e, conse- 


quentemente, temos e 
“DIO Bl) 2A(yplt) | A(O) 
Cop “70 do) “EOP 


Multiplicando Eg. (47) por —[z(t)|>, obtemos a equação diferencial linear fracionária (44). 





(47) 














4.2.1 Equação de Riccati com baixo nível fracionário e coeficientes constantes 


Considere a equação de Riccati (42) com coeficientes constantes e a condição inicial u(0) = 
uo > 0, isto é, 
D “y(t)=ay(t)-by(t)-c, 0<e<1, (48) 


ondea £0,b+0ec 0 são constantes arbitrárias. 








Teorema 4.2 A Eq. (48) tem uma solução da forma y = X, sendo X uma constante se, e somente 
se, X é uma raiz da equação ay? — by —- c = 0. 


O teorema acima permite, entre outras coisas, obter vários exemplos. 
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Exemplo 2 Vamos considerar a equação diferencial fracionária de Riccati, 
Di 'y(e) = Dl) +1, O<e«a, (49) 


com condição inicial y(0) = O. Assim, pelo Teorema (4.2), yp(x) = —1 é uma solução particu- 
lar da Eg. (49). A mudança de variável y(t) = —1 + 7) transforma a Eq. (49) na equação 
diferencial 

Di z(t) +22) =1, (50) 


com z(0) = 1. Aplicando a metodologia acima obtemos a solução geral da Eq. (50) dada por 
1 
zt) = Su +Er-e(—2""8)). (51) 


Finalmente, a solução geral da equação diferencial de Riccati é, assim, dada por 


(1) =-—1+ a 
Po 1+Ere(-28-8) 
AEE) 








— 52 
1+Epe(—21-€) Ro 
Ao aplicar o limite £ — 0" na Eq. (52), recuperamos o resultado clássico 
24 
et —l 
t) = 5—. 53 


5 Conclusão 


A regra da cadeia com baixo nível fracionário justifica alguns versões fracionárias proposto 
por Varalta, Gomes e Camargo (2014) e Tavoni, Mancera e Camargo (2017). Acreditamos 
que nossa abordagem possibilitará a exploração novos caminhos e aplicações usando essa nova 
formulação da regra da cadeia. Além disso, esse método pode ser generalizado para outros tipos 
de Equações Diferenciais Fracionárias, o que será discutido em outros trabalhos. 
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